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Abstract 

Large-scale quantum computation is likely to require massive quantum error correction (QEC). QEC codes and circuits are 
described via the stabilizer formalism, which represents stabilizer states by keeping track of the operators that preserve them. 
Such states are obtained by stabilizer circuits (consisting of CNOT, Hadamard and Phase only) and can be represented compactly 
on conventional computers using fl{n^) bits, where n is the number of qubits |8|. Although techniques for the efficient simulation 
of stabilizer circuits have been studied extensively IJJ, 17J, |8|, techniques for efficient manipulation of stabilizer states are not 
currently available. To this end, we leverage the theoretical insights from 1 1 1 and 1 14] to design new algorithms for: (i) obtaining 
canonical generators for stabilizer states, (n') obtaining canonical stabilizer circuits, and (Hi) computing the inner product between 
stabilizer states. Our inner-product algorithm takes 0{n'^) time in general, but observes quadratic behavior for many practical 
instances relevant to QECC (e.g., GHZ states). We prove that each n-qubit stabilizer state has exactly 4(2" — 1) nearest-neighbor 
stabilizer states, and verify this claim experimentally using our algorithms. We design techniques for representing arbitrary quantum 
states using stabilizer frames and generalize our algorithms to compute the inner product between two such frames. 

I. Introduction 

Gottesman fl) and Knill showed that for certain types of non-trivial quantum circuits known as stabilizer circuits, efficient 
simulation on classical computers is possible. Stabilizer circuits are exclusively composed of stabilizer gates - Controlled-NOT, 
Hadamard and Phase gates (Figure [T^) - followed by measurements in the computational basis. Such circuits are applied to 
a computational basis state (usually |00...0)) and produce output states called stabilizer states. The case of purely unitary 
stabilizer circuits (without measurement gates) is considered often, e.g., by consolidating measurements at the end. Stabilizer 
circuits can be simulated in poly-time by keeping track of a set Pauli operators that stabilizaH^the quantum state. Such stabilizer 
operators uniquely represent a stabilizer state up to an unobservable global phase. Equation^ shows that the number of n-qubit 
stabilizer states grows as 2" therefore, describing a generic stabilizer state requires at least bits. Despite their compact 
representation, stabilizer states can exhibit multi-qubit entanglement and are often encountered in many quantum information 
applications such as Bell states, GHZ states, error-correcting codes and one-way quantum computation. To better understand the 
role stabilizer states play in such applications, researchers have designed techniques to quantify the amount of entanglement [6], 
lfT6l . ifTol in such states and studied relevant properties such as purification schemes ||5], Bell inequalities ||9| and equivalence 
classes fTSl. Efficient algorithms for the manipulation of stabilizer states (e.g., computing the angle between them), can help 
lead to additional insights related to linear-algebraic and geometric properties of stabilizer states. 

In this work, we describe in detail algorithms for the efficient computation of the inner product between stabilizer states. 
We adopt the approach outlined in |r|, which requires the synthesis of a unitary stabilizer circuit that maps a stabilizer state 
to a computational basis state. The work in [IJ shows that, for any unitary stabilizer circuit, there exists an equivalent block- 
structured canonical circuit that applies a block of Hadamard (H) gates only, followed by a block of CNOT (C) only, then 
a block of Phase (P) gates only, and so on in the 7-block sequence H-C-P-C-P-C-H . Using an alternate representation 
for stabilizer states, the work in [14J proves the existence of a (_ff-C-P-CZ)-canonical circuit, where the CZ block consists 
of ControUed-Z (CPHASE) gates. However, no algorithms are known to synthesize such smaller 4-block circuits given an 
arbitrary stabilizer state. In contrast, we describe an algorithm for synthesizing {H-C-C Z-P-H)-cwomca[ circuits given any 
input stabilizer state. We prove that any n-qubit stabilizer state j-i/') has exactly 4(2" — 1) nearest-neighbors - stabilizer states 
|(^) such that KV'I'p)! attains the largest possible value ^ 1. Furthermore, we design techniques for representing arbitrary 
quantum states using stabilizer frames and generalize our algorithms to compute the inner product between two such frames. 

This paper is structured as follows. Section [ll] reviews the stabilizer formalism and relevant algorithms for manipulating 
stabilizer-based representations of quantum states. Section [lll| describes our circuit-synthesis and inner-product algorithms. In 
Section |IV| we evaluate the performance of our algorithms. Our findings related to geometric properties of stabilizer states are 

'An operator U is said to stabilize a state iff U \ip) = \4>)- 
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Fig. 1. (a) Unitary stabilizer gates Hadamard (H), Phase (P) 
and Controlled-NOT (CNOT). 



Fig. 1. (b) The Pauli matrices. 



described in Section |V] In Section VI we discuss stabilizer frames and how they can be used to represent arbitrary states and 
extend our algorithms to compute the inner product between frames. Section [Vll| closes with concluding remarks. 



II. Background and Previous Work 

Gottesman |[8| developed a description for quantum states involving the Heisenberg representation often used by physicists 
to describe atomic phenomena. In this model, one describes quantum states by keeping track of their symmetries rather than 
explicitly maintaining complex vectors. The symmetries are operators for which these states are 1-eigenvectors. Algebraically, 
symmetries form group structures, which can be specified compactly by group generators. It turns out that this approach, also 
known as the stabilizer formalism, can be used to represent an important class of quantum states. 

A. The stabilizer formalism 

A unitary operator U stabilizes a state |-!/)) if \tl>) is a 1-eigenvector of U, i.e., U \'ip) ~ Q, lfT2l . We are interested in 
operators U derived from the Pauli matrices shown in Figure [T]3 The following table lists the one-qubit states stabilized by 
the Pauli matrices. 



(|0)- |1»/V2 
m-^\l))/^/2 
ID 



X : (|0) + |1))/V2 -X 
Y : (|0) +i|l))/\/2 -Y 
Z : |0> -Z \ 

Observe that / stabilizes all states and — / does not stabilize any state. As an example, the entangled state (|00) + |ll))/\/2 
is stabilized by the Pauli operators X®X, —Y®Y, Z®Z and I® I. As shown in Table |I] it turns out that the Pauli matrices 
along with / and the multiplicative factors ±1, ±i, form a closed group under matrix multiplication fTZJ. Formally, the Pauli 
group Qn on n qubits consists of the n-fold tensor product of Pauli matrices, P = i^Pi (g) • • • (8) -Pn such that Pj E {I, X, Y, Z} 
and k e {0, 1, 2, 3}. For brevity, the tensor-product symbol is often omitted so that P is denoted by a string of /, X, Y and Z 
characters or Pauli literals and a separate integer value k for the phase i''. This string-integer pair representation allows us to 
compute the product of Pauli operators without explicitly computing the tensor productsj^] e.g., {— 1 1 X I){iIY 1 1) = —ilYXI. 
Since | Gn \= 4"+^, Qn can have at most logj | Gn \= log2 4"+^ — 2{n + 1) irredundant generators lfT2l . The key idea behind 
the stabilizer formalism is to represent an n-qubit quantum state 1-0) by its stabilizer group S(\ip)) - the subgroup of Qn that 
stabilizes \t/j). As the following theorem shows, if |S'(|7/'))| = 2", the group uniquely specifies 

Theorem 1. For an n-qubit pure state |?/;) and k < n, S(\ip)) = Zj. If k — n, |?/;) is specified uniquely by 5'(|'(/')) and is 
called a stabilizer state. 

Proof (0 To prove that S{\'4>)) is commutative, let P.Q E S{\tp)) such that PQ = If P and Q anticommute, 
-QP = -Q{P lip)) = -Q IV') = - IV') 7^ W- Thus, P and Q cannot both be elements of S{\^)). 

(//) To prove that every element of 5(|?/')) is of degree 2, let P E S{\^p)) such that Pltp) — Observe that P^ = i'/ for 
I E {0, 1, 2, 3}. Since P^ = P{P ^ P IV') = we obtain = 1 and P^ = I. 
(Hi) From group theory, a finite Abelian group with a? = a for every element must be = Z|. 

(iv) We now prove that k < n. First note that each independent generator P E S{\il])) imposes the linear constraint P I?/;) — IV') 
on the 2" -dimensional vector space. The subspace of vectors that satisfy such a constraint has dimension 2""^, or half the space. 
Let gen{\'4))) be the set of generators for S'(|V-')). We add independent generators to gen{\%p)) one by one and impose their 



TABLE I 

Multiplication table for Pauli matrices. Shaded cells 
indicate anticommuting products. 
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linear constraints, to limit to the shared 1-eigenvector. Thus the size of gen{\%lj)) is at most n. In the case |(7en(|-0))| = n, 
the n independent generators reduce the subspace of possible states to dimension one. Thus, is uniquely specified. ■ 
The proof of Theorem [T| shows that S{\'ip)) is specified by only log2 2" = n irredundant stabilizer generators. Therefore, 
an arbitrary n-qubit stabilizer state can be represented by a stabilizer matrix A4 whose rows represent a set of generators 
gi, . ■ . ;.g„ for S{\4>))- (Hence we use the terms generator set and stabilizer matrix interchangeably.) Since each gi is a string 
of n Pauli literals, the size of the matrix is n x n. The phases of each gi are stored separately using a vector of n integers. 
Therefore, the storage cost for is 0(?i^), which is an exponential improvement over the 0(2") cost often encountered in 
vector-based representations. 

Theorem [T] suggests that Pauli literals can be represented using only two bits, e.g., 00 ~ I, 01 = Z, 10 = X and 11 = Y. 
Therefore, a stabilizer matrix can be encoded using an nx2n binary matrix or tableau. The advantage of this approach is that 
this literal-to-bits mapping induces an isomorphism Zj" Gn because vector addition in is equivalent to multiplication 
of Pauli operators up to a global phase. The tableau implementation of the stabihzer formalism is covered in yj, L12J. 

Proposition 2. The number of n-qubit pure stabilizer states is given by 

n-l 



7V(n) = 2"]J(2"-'= + l) 



2(.5+o(l))n2 



(1) 



fe=0 



The proof of Proposition [2] can be found in |1|. An alternate interpretation of Equation [T| is given by the simple recurrence 
relation N{n) — 2(2" + l)N{n — 1) with base case iV(l) = 6. For example, for n — 2 the number of stabilizer states is 60, 
and for n = 3 it is 1080. This recurrence relation stems from the fact that there are 2" + 1 ways of combining the generators 
of N{n — 1) with additional Pauli matrices to form valid n-qubit generators. The factor of 2 accounts for the increase in the 
number of possible sign configurations. Table [II] and Appendix [A| list all two-qubit and three-qubit stabilizer states, respectively. 

Observation 3. Consider a stabilizer state \^p) represented by a set of generators of its stabilizer group S{ \ip)). Recall from 
the proof of Theorem^that, since S{\^p)) = Zj, each generator imposes a linear constraint on \ip). Therefore, the set of 
generators can be viewed as a system of linear equations whose solution yields the 2" basis amplitudes that make up 
Thus, one needs to perform Gaussian elimination to obtain the basis amplitudes from a generator set. 

Canonical stabilizer matrices. Although stabilizer states are uniquely determined by their stabilizer group, the set of generators 
may be selected in different ways. For example, the state \^) — (|00) + |ll))/\/2 is uniquely specified by any of the following 
stabilizer matrices; 



Ml = 
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zz 



M. = 



XX 
-YY 



-YY 

ZZ 



One obtains 7W2 from A^i by left-multiplying the second row by the first. Similarly, one can also obtain A^3 from A^i or 
M.2 via row multiplication. Observe that, multiplying any row by itself yields //, which stabilizes However, // cannot be 



TABLE II 

Sixty two-qubit stabilizer states and their corresponding Pauli generators. Shorthand 

NOTATION REPRESENTS A STABILIZER STATE AS QQ: «!, a2, "3 WHERE di ARE THE NORMALIZED AMPLITUDES 
OF THE BASIS STATES. THE BASIS STATES ARE EMPHASIZED IN BOLD. THE FIRST COLUMN LISTS STATES WHOSE 

GENERATORS DO NOT INCLUDE AN UPFRONT MINUS SIGN, AND OTHER COLUMNS INTRODUCE THE SIGNS, A SIGN 
CHANGE CREATES AN ORTHOGONAL VECTOR, THEREFORE, EACH ROW OF THE TABLE GIVES AN ORTHOGONAL 

BASIS, The CELLS IN DARK GREY INDICATE STABILIZER STATES WITH FOUR NON-ZERO BASIS AMPLITUDES, I,E., 

^ V 4. The Z column indicates the angle between that state and |00), which has 12 

NEAREST-NEIGHBOR STATES (LIGHT GRAY) AND 15 ORTHOGONAL STATES (±), 





State 


Gen'tors 


z 


State 


Gen'tors 


z 


STATE 


Gen'tors 


Z 


State 


Gen'tors 


Z 




1,1,1,1 


IX, XI 


7r/3 


1,-1,1,-1 


-K, XI 


7r/3 


1,1,-1,-1 


IX, 


-XI 


7r/3 


1,-1,-1,1 


-IX, 


-XI 


7r/3 




1,1, j,i 


IX, YI 


7r/3 


1,-1, 


-IX, YI 


7r/3 


1,1,-i, -j 


IX, 


-YI 


7r/3 


1,-1, 


-IX, 


-YI 


7r/3 


_J 


1, 1,0,0 


IX, ZI 


7r/4 


1,-1,0,0 


-IX, ZI 


7r/4 


0,0,1,1 


IX, 


-ZI 


± 


0,0,1,-1 


-IX, 


-ZI 


± 


< 


1,1, l,i 


lY. XI 


ttAS 


1,-i, 1,-i 


-lY, XI 


7r/3 


l,i,-l,-i 


lY, 


-XI 


7r/3 


l,-i,-l,i 


-lY, 


-XI 


7r/3 


< 


lpi(3^j *j — 1 


lY, YI 


7r/3 


1, 1 


-lY, YI 


7r/3 


l,i,-i,l 


lY, 


-YI 


7r/3 


1,-i, -i,-l 


-lY, 


-YI 


7r/3 


w 


l,i,0,0 


lY, ZI 


7r/4 


1, -i,0,0 


-lY, ZI 


7r/4 


0,0, l,i 


lY, 


-ZI 


± 


0,0, 1, -i 


-lY, 


-ZI 




CO 


1,0, 1,0 


IZ, XI 


7r/4 


0, 1,0, 1 


-IZ, XI 


± 


1,0, -1,0 


IZ, 


-XI 


7r/4 


0,1,0,-1 


-IZ, 


-XI 






1,0, i,0 


IZ, YI 


7r/4 


0, l,0,i 


-IZ, YI 




1,0, -i,0 


IZ, 


-YI 


7r/4 


0, 1,0, -i 


-IZ, 


-YI 






1,0,0,0 


IZ, ZI 





0,1,0,0 


-IZ, ZI 




0,0,1,0 


IZ, 


-ZI 


L 


0,0,0,1 


-IZ, 


-ZI 


L 


Q 


0, 1, 1,0 


XX, YY 


L 


1,0,0,-1 


-XX, YY 


7r/4 

± 


1,0,0,1 


XX, 


-YY 


7r/4 


0,1,-1,0 


-XX, 


-YY 


L 


_j 


1,0, 0,i 


XY, YX 


7I-/4 


0,1,1,0 


-XY, YX 


0,1, -i,0 


XY, 


-YX 


± 


1,0,0, -i 


-XY, 


-YX 


7r/4 


o 
z 


1,1,1,-1 


XZ, ZX 


7r/3 


1,1,-1,1 


-XZ, ZX 


7r/3 


1,-1,1,1 


XZ, 


-ZX 


7r/3 


1,-1,-1,-1 


-XZ, 


-ZX 


7r/3 


< 


1,1, 1,-i 


XZ, ZY 


7r/3 


1,1, -l,i 


-XZ, ZY 


7r/3 


l,-i,l,i 


XZ, 


-ZY 


7r/3 


l,-i,-l,-j 


-XZ, 


-ZY 


7r/3 


z 


l,l,i,-i 


YZ, ZX 


7r/3 


1, 1, 


-YZ, ZX 


7r/3 


1,-1, i,z 


YZ, 


-ZX 


7r/3 


1,-1,-i, -j 


-YZ, 


-ZX 


7r/3 


w 


l,i,i,l 


YZ, ZY 


7r/3 


l,i,-i,-l 


-YZ, ZY 


7r/3 


l,-i,i,-l 


YZ, 


-ZY 


7r/3 


l,-i,-i, 1 


-YZ, 


-ZY 


7r/3 



3 



LX, Y, Z literals 



'Son 



I, Z literals 



I, Z literals 



/ literals only 



Fig. 2. Canonical (row-reduced echelon) form for stabilizer matrices. 
The X-block contains a minimal set of rows with XIY literals. The 
rows with Z literals only appear in the Z-block. Each block is arranged 
so that the leading non-/ literal of each row is strictly to the right of 
the leading non-/ literal in the row above. The number of Pauli (non-/) 
literals in each block is minimal. 



used as a stabilizer generator because it is redundant and carries no information about the structure of \'ip). This also holds true 
in general for of any size. Any stabilizer matrix can be rearranged by applying sequences of elementary row operations in 
order to obtain a particular matrix structure. Such operations do not modify the stabilizer state. The elementary row operations 
that can be performed on a stabilizer matrix are transposition, which swaps two rows of the matrix, and multiplication, which 
left-multiplies one row with another. Such operations allow one to rearrange the stabilizer matrix in a series of steps that 
resemble Gauss-Jordan elimination]^ Given an n x n stabilizer matrix, row transpositions are performed in constant tim^ 
while row multiplications require 8(n) time. Algorithm [T] rearranges a stabilizer matrix into a row-reduced echelon form that 
contains: (/) a minimal set of generators with X and Y literals appearing at the top, and (//) generators containing a minimal set 
of Z literals only appearing at the bottom of the matrix. This particular stabilizer-matrix structure, shown in Figure [2] defines a 
canonical representation for stabilizer states. The algorithm iteratively determines which row operations to apply based on the 
Pauli (non-/) literals contained in the first row and column of an increasingly smaller submatrix of the full stabilizer matrix. 
Initially, the submatrix considered is the full stabilizer matrix. After the proper row operations are applied, the dimensions of 
the submatrix decrease by one until the size of the submatrix reaches one. The algorithm performs this process twice, once 
to position the rows with X{Y) literals at the top, and then again to position the remaining rows containing Z literals only at 
the bottom. Let i e {1, . . . , n} and j e {1, . . . , n} be the index of the first row and first column, respectively, of submatrix 
A. The steps that achieve the upper-triangular of the row-echelon form shown in Figure |2] are as follows. 

1. Let fc be a row in A whose j*'* literal is X{Y). Swap rows k and i such that k is the 1** row of A. Decrease the height 
of A by one (i.e., increase i). 

2. For each row m G {0, . . . , n}, m ^ i that has an X(Y) in column j, use row multiplication to set the j*^ literal in row 
m to / or Z. 

3. Decrease the width of A by one (i.e., increase j). 

To construct the lower-triangular of the matrix, one simply executes the same process with the following differences: (/) 
step 1 looks for rows that have a Z or Y literals in column j, and (//) the height i of A ranges from {r + 1, . . . n}, where r is 
the number of rows that have X oi Y literals. Since Algorithm [l] looks at all entries in the matrix and performs a constant 
number of row multiplications each time, the runtime of the algorithm is 0{n^). 

Stabilizer circuit simulation. The computational basis states are stabilizer states that can be represented using the following 
stabilizer-matrix structure. 

Definition 4. A stabilizer matrix is in basis form if it has the following structure. 



± 
± 

± 



Z I ■■■ I 
I Z ■■■ I 

I I ■■■ z 



In this matrix form, the ± sign of each row along with its corresponding -literal designates whether the state of the j*-^ 
qubit is |0) (+) or |1) (— ). Suppose we want to simulate circuit C. Stabilizer-based simulation first initializes to specify 
some basis state To simulate the action of each gate [/ e C, we conjugate each row gi of by C/j^ We require that 
UgiW maps to another string of Pauli literals so that the resulting stabilizer matrix Ai' is well-formed. It turns out that the 
Hadamard, Phase and CNOT gates (Figure [T^) have such mappings, i.e., these gates conjugate the Pauli group onto itself ||8], 
lfT2l . Table III lists the mapping for each of these gates. 



Since Gaussian elimination essentially inverts the n X 2n matrix, this could be sped up to 0(n^*^^^) time by using fast matrix inversion algorithms. However, 0{n'^)- time 
Gaussian elimination seems more practical. 

^Storing pointers to rows facilitates 0(l)-time row transpositions - one simply swaps relevant pointers. 

^Since gi \ip) = \ip), the resulting state U \ip) is stabilized by UgiU^ because (UgiU^)U \ip) = Ugi \tp) = U lii). 
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Algorithm 1 Canonical form reduction for stabilizer matrices 



Input: Stabilizer matrix A4 for S(\ip)) with rows Ri, . . . , _R„ 
Output: A4 is reduced to row-echelon form 

=> ROVSVAP{M, i,j) swaps rows Ri and Rj of M 

=> R0WMULT(A4, j, j) left-multiplies rows i?i and Rj, returns updated Ri 



i ^ 1 

for j G {1, . . . , n} do 

k index of row Rks{i....,n} with j*'' literal set to X(y) 
if k exists then 

R0WSWAP(>1,i,fc) 

for me {0, . . . , n} do 

if j"* literal of i?™ is X or F and m ^ i then 

i?„ = ROWMULT(X, Ri,Rm) 
end if 
end for 
i i + 1 
end if 
end for 

for j G {1, . . . , n} do 

k index of row Rks{i,...,n} with j"' literal set to Z 
if k exists then 

R0WSWAP(>1,i,fc) 
for m e {0, . . . , n} do 

if J*'' literal of Rm. is Z or y and m ^ i then 

7?„ = R0WMULT(>1, J?i, Rm) 
end if 
end for 
i i + 1 
end if 
end for 



> Setup X block 



> Gauss-Jordan elimination step 



> Setup Z block 



> Gauss-Jordan elimination step 



For example, suppose we simulate a CNOT operation on lip) = (|00) + |ll))/v2 using A4, 



^= zz 



XI 
IZ 



C NOT J — 

One can verify that the rows of M.' stabilize j?/') > (|00) + |10))/v2 as required. 

Since Hadamard, Phase and CNOT gates are directly simulated using stabilizers, these gates are commonly called stabilizer 
gates. They are also called Clijford gates because they generate the Clifford group of unitary operators. We use these names 
interchangeably. Any circuit composed exclusively of stabilizer gates is called a unitary stabilizer circuit. Table [Hi] shows that 
at most two columns of are updated when one simulates a stabilizer gate. Thus, such gates are simulated in <d{n) time. 

Theorem 5. An n-qubit stabilizer state can be obtained by applying a stabilizer circuit to the jO)**" basis state. 

Proof: The work in H] represents the generators using a tableau, and then shows how to construct a unitary stabihzer 
circuit from the tableau. We refer the reader to HI Theorem 8] for details of the proof. ■ 

Corollary 6. An n-qubit stabilizer state j^/j) can be transformed by stabilizer gates into the |0)^" basis state. 

Proof: Since every stabilizer state can be produced by applying some unitary stabilizer circuit C to the 10)*^" state, it 
suffices to reverse C to perform the inverse transformation. To reverse a stabilizer circuit, reverse the order of gates and replace 
every P gate with PPP. ■ 

The stabilizer formalism also admits one-qubit measurements in the computational basis |8|. However, the updates to 
for such gates are not as efficient as for stabilizer gates. Note that any qubit j in a stabilizer state is either in a |0) (|1)) state 
or in an unbiasecj^ superposition of both. The former case is called a deterministic outcome and the latter a random outcome. 
We can tell these cases apart in Q{n) time by searching for X or Y literals in the j*'' column of Ai. If such literals are 
found, the qubit must be in a superposition and the outcome is random with equal probability {p(0) = p(l) = .5); otherwise 
the outcome is deterministic (p(0) = 1 or p{l) — 1). 

^An arbitrary state with computational basis decomposition Y2k=o l'^) ^^^^ ^ unbiased if for all ^ and Xj ^0, |Aip — |Aj p. Otherwise, the state is 
biased. One can verify that none of the stabilizer gates produce biased states. 
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TABLE III 

Conjugation of the Pauli-group elements by the stabilizer gates (12). 
For the CNOT case, subscript 1 indicates the control and 2 the target. 



Gate 


Input 


Output 




X 


Z 


H 


Y 


-Y 




Z 


X 




X 


Y 


P 


Y 


-X 




Z 


Z 



Gate 


Input 


Output 


CNOT 


I1X2 
X1I2 
I1Y2 
Y1I2 

hZ2 


hX2 

X\ X2 
Z1Y2 
Y1X2 
Z1Z2 
Zih 



Random case: one flips an unbiased coin to decide the outcome and then updates to make it consistent with the outcome 
obtained. This requires at most n row muhiplications leading to 0{n?) runtime HI, fl2|. 

Deterministic case: no updates to are necessary but we need to figure out whether the state of the qubit is |0) or |1), i.e., 
whether the qubit is stabiHzed by Z or -Z, respectively. One approach is to apply Algorithm [T] to put M. in row-echelon form. 
This removes redundant literals from M in order to identify the row containing & Z in its j position and / everywhere else. 
The ± phase of this row decides the outcome of the measurement. Since this approach is a form of Gaussian elimination, it 
takes 0{ii?) time in practice. 

Aaronson and Gottesman 1 1 1 improved the runtime of deterministic measurements by doubling the size of M to include n 
destabilizer generators in addition to the n stabilizer generators. Such destabilizer generators help identify exactly which row 
multiplications to compute in order to decide the measurement outcome. This approach avoids Gaussian elimination and thus 
deterministic measurements are computed in 0{v?) time. 

III. Inner-product and circuit-synthesis algorithms 
Given (V'lv') = re*", we normalize the global phase of j?/') to ensure, without loss of generality, that G M+. 

Theorem 7. Let S'(|'0)) and S{\lp)) be the stabilizer groups for |?/;) and \ip), respectively. If there exist P £ S'd?/')) and 
Q G S{\(p)) such that P = -Q, then _L \ip). 

Proof Since is a 1-eigenvector of P and \(p) is a (—1) -eigenvector of P, they must be orthogonal. ■ 

Theorem 8. [1| Let 1-0), \ip) be non-orthogonal stabilizer states. Let s be the minimum, over all sets of generators {Pi, . . . , P„} 
for S{\ij)) and {Qi, ■ ■ ■ ,Qn} fof S{\ip)), of the number of dijferent i values for which P^ 7^ Qi. Then, = 2^^*/^. 

Proof: Since {ipltp) is not affected by unitary transformations U, we choose a stabilizer circuit such that U — \b), where 
1 6) is a basis state. For this state, select the stabilizer generators Ai of the form / . . . IZI . . . I. Perform Gaussian elimination 
on Ai to minimize the incidence of Pi 7^ Qi. Consider two cases. If U \ (p) \b) and its generators contain only I/Z literals, 
then U \ tp) _L U {tp), which contradicts the assumption that and \ip) are non-orthogonal. Otherwise, each generator of U \ (p) 
containing X/Y literals contributes a factor of 1/ \/2 to the inner product. ■ 

Synthesizing canonical circuits. A crucial step in the proof of Theorem |8] is the computation of a stabilizer circuit that brings 
an n-qubit stabilizer state to a computational basis state \b). Consider a stabilizer matrix Ai that uniquely identifies 
A4 is reduced to basis form (Definition |4| by applying a series of elementary row and column operations. Recall that row 
operations (transposition and multiplication) do not modify the state, but column (Clifford) operations do. Thus, the column 
operations involved in the reduction process constitute a unitary stabilizer circuit C such that C\4>) = \b), where \b) is a basis 
state. Algorithm |2] reduces an input stabilizer matrix Ai to basis form and returns the circuit C that performs such a mapping. 

Definition 9. Given a finite sequence of quantum gates, a circuit template describes a segmentation of the circuit into blocks 
where each block uses only one gate type. The blocks must correspond to the sequence and be concatenated in that order. For 
example, a circuit satisfying the H-C-P template starts with a block of Hadamard {H) gates, followed by a block of CNOT 
(C) gates, followed by a block of Phase (P) gates. 

Definition 10. A circuit with a template structure consisting entirely of CNOT, Hadamard and Phase blocks is called a 
canonical stabilizer circuit. 

Canonical forms are useful for synthesizing stabilizer circuits that minimize the number of gates and qubits required to 
produce a particular computation. This is particularly important in the context of quantum fault-tolerant architectures that are 
based on stabilizer codes. Given any stabilizer matrix. Algorithm |2] synthesizes a 5-block canonical circuit with template H- 
C-CZ-P-H (Figure |3]-a), where the CZ block consists of Controlled-Z (CPHASE) gates. Such gates are stabilizer gates since 
CPHASEi j =HjCNOTi jHj (Figure |3]-b). In our implementation, such gates are simulated directly on the stabilizer The work 
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Fig. 3. (a) Template structure for the basis-normalization circuit synthesized by Algorithm|2] The input is an arbitrary 
stabilizer state while the output is a basis state \bi, . . . , b„), where bi, . . . ,b„ G {0, 1}^. (b) ControUed-Z gates 
used in the CPHASE block. CPHASE gates can be implemented directly or using the equivalence shown here. 

in Q establishes a longer 7-hloc^H-C-P-C-P-C-H canonical-circuit template. The existence of a H-C-P-C Z template is 
proven in fl4l but no algorithms are known for obtaining such 4-block canonical circuits given an arbitrary stabilizer state. 

We now describe the main steps in Algorithm |2] For simplicity, the updates to the phase array under row and column 
operations will be left out of our discussion as such updates do not affect the overall execution of the algorithm. 

1. Reduce M to canonical form. 

2. Use row transposition to diagonalize M. For j E {1, . . . , n}, if the diagonal literal Mjj = Z and there are other Pauli 
(non-J) literals in the row (qubit is entangled), conjugate by Hj. Elements below the diagonal are Z/I literals. 

3. For each above-diagonal element Aij^k = XIY, conjugate by CNOTj ^. Elements above the diagonal are now IIZ literals. 

4. For each above-diagonal element Mj_k — Z, conjugate by CPHASEj ^. Elements above the diagonal are now / literals. 

5. For each diagonal literal A^jj = Y, conjugate by P^. 

6. For each diagonal literal — X, conjugate by Hj. 

7. Use row multiplication to eliminate trailing Z literals below the diagonal and arrive at basis form. 

Proposition 11. For an n x n stabilizer matrix M, the number of gates in the circuit C returned by Algorithm^is 0{n^). 

Proof: The number of gates in C is dominated by the CNOT and CPHASE blocks, which have 0{n?) gates each. This 
agrees with previous results regarding the number of gates needed for an n-qubit stabilizer circuit in the worst case lIS), fH. ■ 

Observe that, for each gate added to C, the corresponding column operation is applied to A4. Therefore, since column 
operations run in Q{n) time, it follows from Proposition 11 that the runtime of Algorithm |2] is 0{n^). 

Canonical stabilizer circuits that follow the 7-block template structure from ||T| can be optimized to obtain a tighter bound 
on the number of gates. As in our approach, such circuits are dominated by the size of the CNOT blocks, which contain 0{n^) 
gates. The work in |13 | shows that that any CNOT circuit has an equivalent CNOT circuit with only 0{n^ / \ogn) gates. Thus, 
one simply applies such techniques to each of the CNOT blocks in the canonical circuit. It is an open problem whether one 
can apply the techniques from llT3l directly to CPHASE blocks, which would facilitate similar optimizations to our proposed 
5-block canonical form. 

Inner-product algorithm. Let and \(t>) be two stabilizer states represented by stabilizer matrices Ai"^ and Ai'^, respectively. 
Our approach for computing the inner product between these two states is shown in Algorithm[3] Following the proof of Theorem 
[s] Algorithm [2] is applied to A4'^ in order to reduce it to basis form. The stabilizer circuit generated by Algorithm [2] is then 
applied to A^^in order to preserve the inner product. Then, we minimize the number of X and Y literals in A^"^ by applying 
Algorithm [1] Finally, each generator in A^"^ that anticommutes with Ai '^' (since Al"^ is in basis form, we only need to check 
which generators in Ai'^ have X or Y literals) contributes a factor of l/\/2 to the inner product. If a generator in Ai'^, say 
Qi, commutes with Ai"^, then we check orthogonality by determining whether Qi is in the stabilizer group generated by Ai"^'. 
This is accomplished by multiplying the appropriate generators in AA"^ such that we create Pauli operator R, which has the 
same literals as Qi, and check whether R has an opposite sign to Qi. If this is the case, then, by Theorem [7] the states are 
orthogonal. Clearly, the most time-consuming step of Algorithm |3] is the call to Algorithm |2] therefore, the overall runtime 



is 0{n ). However, as we show in Section IV the performance of our algorithm depends strongly on the stabilizer matrices 



considered and exhibits quadratic behaviour for certain stabilizer states. 

IV. Empirical Validation 

We implemented our algorithms in C++ and designed a benchmark set to validate the performance of our inner-product 
algorithm. Recall that the runtime of Algorithm [2] is dominated by the two nested for-loops (lines 20-35). The number of times 

^Theorem 8 in 1 1 1 actually describes an 11-step canonical procedure. However, the last four steps pertain to reducing destabilizer rows, which we do not 
consider in our approach. 



7 



Algorithm 2 Synthesis of basis normalization circuit 
Input: Stabilizer matrix A4 for S{\i})) with rows , . . . , Rn 

Output: (0 Unitary stabilizer circuit C such that C lip) equals basis state \b), and [ii) reduce A4 to basis form 
=> GAUSS reduces M to canonical form (Figure |2]( 
=> RO\lSVAP{M,i,j) swaps rows Ri and Rj of M 
=> R0WMULT(A1, j, j) left-multiphes rows Ri and _Rj, returns updated Ri 
C0NJ(A1,aj) conjugates j*'' column of M by Clifford sequence a 



GAUSS > Set M to canonical form 

C ^ 
i ^ 1 

for j £ {1, . . . , n} do l> Apply block of Hadamard gates 

k index of row Rke{i,...,n} with j*^ literal set to X or y 
if k exists then 

ROVSVkP{M,i,k) 
else 

k2 index of last row -Rfc2e{', ■■■,«} ^i^^ j"" literal set to Z 
if fc2 exists then 

R0WSWAP(X,i,fc2) 

if Ri has X, Y or Z literals in columns {j + 1, . . . ,n} then 
C0NJ(A1,Hj) 
C C U Hj 
end if 
end if 
end if 
i i + 1 
end for 

for j G {1, . . . , n} do > Apply block of CNOT gates 

for fc G {j + 1, . . . , n} do 

if fc*'^ literal of row Rj is set to A or y then 
C0NJ(A1,CNOTj,fe) 
C^CU CNOTj,fe 
end if 
end for 
end for 

for j £ {1, . . . , n} do > Apply a block of ControUed-Z gates (Figure [sj)) 

for fc e {j + 1, . . . , n} do 

if fc*** literal of row Rj is set to Z then 
C0NJ(>1,CPHASEj,fc) 
C ^ C U CPHASEj,fc 
end if 
end for 
end for 

for j G {1, . . . , n} do l> Apply block of Phase gates 

if j"' literal of row Rj is set to Y then 

CONJ{M,Pj) 

c^cuPj 

end if 
end for 

for j G {1, . . . , n} do l> Apply block of Hadamard gates 

if j*'' literal of row Rj is set to X then 
C0NJ(7V(,Hj) 
C^CUHj 
end if 
end for 

for j G {1, . . . , n} do > Eliminate trailing Z literals to ensure basis form (Definition [4j 

for fc G {j + 1, . . . , n} do 

if j"' literal of row Rt is set to Z then 

Rk= RDWMULT {M,Rj,Rk) 
end if 
end for 
end for 
return C 
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Algorithm 3 Inner product for stabilizer states 

Input: Stabilizer matrices (i) M"^ for \ip) with rows Pi, ... , P„, and (ii) M"^ for |(^) witfi rows Qi, . . . , Q„ 
Output: Inner product between and 

=> BASISNORMCIRC(AI) reduces M to basis form, i.e, C \ip) = |fe), where \h) is a basis state, and returns C 

=> C0NJ(7\4,C) conjugates 7\4 by Clifford circuit C 

=> GAUSS reduces M to canonical form (Figure [2]l 

LEFTMULT(P, Q) left-multiplies Pauli operators P and Q, and returns the updated Q 

C ^ BASISN0RMCIRC(>1'*) 
CONJ(X'^,C) 
GAUSS (yvc*) 
fc ^ 

for each row Qi G A^'* do 

if Qi contains X oi Y literals then 

fc ^ fc + 1 

else 

P^ /®" 

for each Z literal in Qi found at position j do 

R -h- LEFTMULT(Pj,P) 
end for 

if P = — Qi then 

return 
end if 
end if 
end for 
return 2 '^''^ 



> Apply Algorithm [2] to M"^ 
t> Compute C \ (j>) 
t> Set A^** to canonical form 



t> Check orthogonahty, i.e., Qi ^ S(\b)). 

> By Theorem |7] 
t> By Theorem [s] 



these loops execute depends on the amount of entanglement in the input stabilizer state. In turn, the number of entangled 
qubits depends on the the number of CNOT gates in the circuit C used to generate the stabilizer state C |0®") (Theorem [sjl. 
By a simple heuristic argument 1 1 1, one generates highly entangled stabilizer states as long as the number of CNOT gates in C 
is proportional to nlgn. Therefore, we generated random n-qubit stabilizer circuits for n S {20,40, . . . , 500} as follows: fix 
a parameter /3 > 0; then choose /3\n\og2 n] unitary gates (CNOT, Phase or Hadamard) each with probability 1/3. Then, each 
random C is applied to the |00 . . . 0) basis state to generate random stabilizer matrices (states). The use of randomly generated 
benchmarks is justified for our experiments because (/) our algorithms are not explicitly sensitive to circuit topology and (//) 
random stabilizer circuits are considered representative [ 1 1 1 . For each n, we applied Algorithm |3] to pairs of random stabilizer 
matrices and measured the number of seconds needed to compute the inner product. The entire procedure was repeated for 
increasing degrees of entanglement by ranging (3 from 0.6 to 1.2 in increments of 0.1. Our results are shown in Figure |4]-a. 

The runtime of Algorithm [3] appears to grow quadratically in n when (3 = 0.6. However, when the number of unitary gates 
is doubled (/3 = 1.2), the runtime exhibits cubic growth. Therefore, Figure |4]-a shows that the performance of Algorithm [3] 
is highly dependent on the degree of entanglement in the input stabilizer states. Figure Hl-b shows the average size of the 
basis-normalization circuit returned by the calls to Algorithm |2] As expected (Proposition|lT|, the size of the circuit grows 
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Fig. 4. Average runtime for Algorithm|3]to compute the inner product between two random n-qubit stabilizer states. 
The stabilizer matrices that represent the input states are generated by applying /3nlog2 n unitary gates to |0®"). 
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Fig. 5. Average runtime for Algoritlim 3 to compute the inner product between (a) 
state I (f>) and (b) the n-qubit GHZ state ana random stabilizer state \<f>). 



and random stabilizer 



quadratically in n. Figure |5] shows the average runtime for Algorithm |3] to compute the inner product between: (/) the all- 
zeros basis state and random n-qubit stabilizer states, and (//) the n-qubit GHZ stat^ and random stabilizer states. GHZ 
states are maximally entangled states that have been realized experimentally using several quantum technologies and are often 
encountered in practical applications such as error-correcting codes and fault-tolerant architectures. Figure [5] shows that, for 
such practical instances. Algorithm |3] can compute the inner product in roughly 0{n^) time (e.g. {GHZ\Q)). However, without 
apriori information about the input stabilizer matrices, one can only say that the performance of Algorithm |3] will be somewhere 
between quadratic and cubic in n. 

V. Nearest-neighbor stabilizer states 

We used Algorithm [3] to compute the inner product between |00) and all two-qubit stabiUzer states. Our results are shown 
in Table |ll] We leveraged these results to formulate the following properties related to the geometry of stabilizer states. 



Definition 12. Given an arbitrary state with 
attains the largest possible value ^ 1. 

Proposition 13. Consider two orthogonal stabilizer states \a) and 
state. Then \^) is a nearest stabilizer state to \a) and \j3). 

Proof: Since stabilizer states are unbiased, KV'lo:) 
Thus, IV') is a nearest stabilizer state to \a) and |/3). 



1, a stabilizer state \tf) is a nearest stabilizer state to if K'f/'li/?)! 

whose unbiased superposition is also a stabilizer 
KV'I/?)! = By Theorem jsj this is the largest possible value ^ 1. 



Lemma 14. For any two stabilizer states, the numbers of nearest-neighbor stabilizer states are equal. 

Proof: By Corollary |6] any stabilizer state can be mapped to another stabilizer state by a stabilizer circuit. Since the 
operators effected by these circuits are unitary, inner products are preserved. ■ 

Lemma 15. Let \^) and \lp) be orthogonal stabilizer states such that \ip) — P\'ip) where P is an element of the Pauli group. 
Then ^'^'y^'^^ is a stabilizer state. 

Proof: Suppose \il>) = {gk)k=i,2,. 



is generated by elements of the ri-qubit Pauli group. Let 



if[P,.9fc]=0 

1 otherwise 

and write \ip) = {{— 1)^ gk). Conjugating each generator gk by P we see that \ip) is stabilized by gk). Let 

(respectively Xk) denote the Pauli operator Z (X) acting on the fc*'' qubit. By Corollary |6] there exists an element L of the 
n-qubit Clifford group such that L\ij) = |0)'^" and L\ip) = {LPL^)L\'4j) = i*|/(l)/(2) . . . /(n)). The second equaUty follows 
from the fact that LPL^ is an element of the Pauli group and can therefore be written as i*X{v)Z{u) for some t E {0, 1, 2, 3} 

and Therefore, + Lt(|0)^" + |/(l)/(2) . . . /(n))) 



^/2 



V2 



An n-qubit GHZ state is an equal superposition of the all-zeros and all-ones states, i.e. 



|0»") + |1«") 
V2 
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The state in parenthesis on the right-hand side is the product of an all-zeros state and a GHZ state. Therefore, the sum is 
stabilized by S' = {Szeroi Sghz)L where Szero = {Zi,i G {k\f{k) — 0}) and Sghz is supported on {k\f{k) = 1} and 
equals ((-l)*/2XX...X,Vi if t = mod 2 or ((-l)(*-i)/2yr . . . F, Vi if i = 1 mod 2. ■ 

Theorem 16. For any n-qubit stabilizer state \^), there are 4(2" — 1) nearest-neighbor stabilizer states, and these states can 



be produced as described in Lemma 15. 



Proof: The all-zeros basis amplitude of any stabilizer state \^) that is a nearest neighbor to 10)*^" must be oc l/\/2. 
Therefore, ji/;) is an unbiased superposition of lO)**" and one of the other 2" — 1 basis states, i.e., = — — ^J^^^ — , where 



P €Qn such that P |0)®" 7^ a |0)^". As in the proof of Lemma 15 we have = '"^ v^''''^^ where \ip) is a basis state 
and t G {0, 1, 2, 3}. Thus, there are 4 possible unbiased superpositions, and a total of 4(2" — 1) nearest stab 
|0)^" is a stabilizer state, all stabilizer states have the same number of nearest stabilizer states by Lemma 

Table [n] shows that |00) has 12 nearest-neighbor states. We computed inner products between all-pairs of 2-qubit stabilizer 
states and confirmed that each had exactly 12 nearest neighbors. We used the same procedure to verify that all 3-qubit stabilizer 
states have 28 nearest neighbors. We verified the correctness of our algorithm by comparing against inner product computations 
based on explicit basis amplitudes. 

VI. Stabilizer frames 

Given an n-qubit stabilizer state there exists an orthonormal basis including and consisting entirely of stabilizer 
states. Using Theorem [7] one can generate such a basis from the stabilizer representation of \^]). Observe that, one can create 
a state \lp) that is orthogonal to by changing the signs of an arbitrary non-empty subset of generators of S{\'4))), i.e., by 
permuting the phase vector of the stabilizer matrix for \'4>). Moreover, selecting two different subsets will produce two mutually 
orthogonal states. Thus, one can produce 2" — 1 additional orthogonal stabilizer states. Such states, together with !?/>), form 
an orthonormal basis. This is illustrated by Table |ll] were each row constitutes an orthonormal basis. 

Definition 17. A stabilizer frame is a set of /c < 2" stabilizer states that forms an orthonormal basis {1-01} , . . . , 1^/;^)} and 
spans a subspace of the n-qubit Hilbert space. F is represented by a pair consisting of (/) a stabilizer matrix and (//) a set 
of k distinct phase vectors (Tj{A4),j G {1, . . . , k}. The size of the frame, which we denote by is equal to k. 

Stabilizer frames are useful for representing arbitrary quantum states and for simulating the action of stabilizer circuits on 
such states. Let a. — {ai, . . . ,ak) G C*^ be the decomposition of the arbitrary n-qubit state |(/)) onto the basis {Itpi) , . . . , \ipk)} 
defined by i.e., \(f>) = X]^=i "^fc li^i}- Furthermore, let [/ be a stabilizer gate. To simulate C/|(/)), one simply rotates the 
basis defined by J- to get the new basis {U j-fAi) , ■ ■ ■ ,U \4'k)}- This is accomplished with the following two-step process: (/) 
update the stabilizer matrix A4 associated with as per Section |11-A[ (//) iterate over the phase vectors in and update each 
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accordingly (Table llli. The second step is linear in the number of phase vectors as only a constant number of elements in 
each vector needs to be updated. Also, ct may need to be updated, which requires the computation of the global phase of 
each U lipi)- Since the stabilizer does not maintain global phases directly, each a,; can be updated as follows: (/) reduce M to 
canonical form, (;/) obtain a basis state \b) from A4 and store its non-zero amplitude /3, (///) obtain U \b) from UMU^ and 
store its non-zero amplitude 7, and (/v) compute — ai * Observe that, all the above processes take time polynomial 
in k, therefore, if fc = poly{n), U \(j>) can be simulated efficiently on a classical computer via frame-based simulation. 



Inner product between frames. We now discuss how to use our algorithms to compute the inner product between arbitrary 
quantum states. Let |0) and \(p) be quantum states represented by the pairs < J^"^, a = (ai, . . . ,0;^) > and < J^'^, (3 = 
(/3i, . . . , (3i) >, respectively. The following steps compute 

1. Apply Algorithm [2] to A4'^ (the stabilizer matrix associated with J^'^) to obtain basis-normaUzation circuit C. 

2. Rotate frames T'^ and J^"^ by C as outlined in our previous discussion. 

3. Reduce A^"^ to canonical form (Algorithm [T} and record the row operations applied. Apply the same row operations to 
each phase vector af , i G {1, . . . , A:} in T'^. Repeat this step for Ai'-^ and the phase vectors in T'^. 

4. Let Aif denote that the leading-phases of the rows in A^"^ are set equal to af. Similarly, AiJ denotes that the phases 
of Ai"^ are equal to crj. Furthermore, let 5{Aif , Ai'j) be the function that returns if the orthogonality check from 
Algorithm |3] (lines 9-15) returns 0, and 1 otherwise. The inner product is computed as, 

k I 

where s is the number of rows in Aiu, that contain X 01 Y literals. 



11 



Prior work on representation of arbitrary states using the stabilizer formalism can be found in |[T|. The authors propose an 
approach that represents a quantum state as a sum of density-matrix terms. Our frame-based technique offers more compact 
storage < 2" whereas a density matrix may have 4" non-zero entries) but requires more sophisticated book-keeping. 

VII. Conclusion 

The stabilizer formalism facilitates compact representation of stabilizer states and efficient simulation of stabilizer circuits. 
Stabilizer states appear in many different quantum-information applications, and their efficient manipulation via geometric and 
linear-algebraic operations may lead to additional insights. To this end, we study algorithms to efficiently compute the inner 
product between stabilizer states. A crucial step of this computation is the synthesis of a canonical circuit that transforms a 
stabilizer state into a computational basis state. We designed an algorithm to synthesize such circuits using a 5-block template 
structure and showed that these circuits contain 0{n^) stabilizer gates. We analysed the performance of our inner-product 
algorithm and showed that, although its runtime is 0{n^), there are practical instances in which it runs in linear or quadratic 
time. Furthermore, we proved that an n-qubit stabilizer state has exactly 4(2" — 1) nearest-neighbor states and verified this 
result experimentally. Finally, we designed techniques for representing arbitrary quantum states using stabilizer frames and 
generalize our algorithms to compute the inner product between two such frames. 
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Appendix A 

The 1080 three-qubit stabilizer states 

Shorthand notation represents a stabihzer state as ^i, 0^2, oiz where aj are the normalized amphtudes of the basis states. 
Basis states are emphasized in bold. The Z column indicates the angle between that state and |000), which has 28 nearest- 
neighbor states and 315 orthogonal states (_L). 
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